Optimal tuning of solid-state quantum gates: A universal two-qubit gate 
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We present a general route to reduce inhomogeneous broadening in nanodevices due to 1// noise. 
We apply this method to a universal two-qubit gate and demonstrate that for selected optimal 
couplings, a high-efficient gate can be implemented even in the presence of 1 // noise. Entanglement 
degradation due to interplay of 1// and quantum noise is quantified via the concurrence. A charge- 
phase vi — SWAP gate for spectra extrapolated from single qubit experiments is analyzed. 
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Intense research on solid state nanodevices during the 
last decade lead to observation of fundamental quantum 
phenomena at the nanoscale. Combining quantum coher- 
ence with the existing intcgrated-circuit fabrication tech- 
nology makes nanocircuits very promising for quantum 
computing. In particular, a variety of high-fidelity single 
qubit gates on a superconducting platform are nowadays 
available^— . Controlled generation of entangled states 
and preserving quantum correlations represent a timely 
and critical issue. Identification of strategies to counter- 
act physical processes detrimental to quantum coherent 
behavior is a fundamental step towards this goal. 

Noise with 1// spectrum is ubiquitous in nanodevices 
and represents a serious limitation for their coherence 
properties^—. Because of low frequency fluctuations of 
the device eigenenergies, the average (implied by quan- 
tum measurement) of signals occurring at slightly differ- 
ent frequencies, is defocuscd. This effect, analogous to 
inhomogeneous broadening^, is a signature of 1/f noise, 
whatever its microscopic origin. It results in a typical 
algebraic decay of coherent oscillations^, as observed in 
single-qubit gates of various design and materials 2 -^. 

One way to defeat inhomogeneities is exploiting NMR 
decoupling procedures^. For instance, echo-type proto- 
cols may considerably reduce defocusing in single qubit 
rotations 2 ^^, but they also limit external control on 
quantum state processing. Moreover, extension to multi- 
qubit gates may require a large fraction of quantum op- 
erations devoted to effective decoupling procedures, lim- 
iting scalability. A simpler but effective strategy is to 
operate the device at "optimal points" , characterized by 
minimal sensitivity of the relevant splittings to variations 
of the control parameters. Partial reduction of defocusing 
of single-qubit operations by operating at such working 
point has been successfully demonstrate d 1 ' 2 ' 10 . 

In this Brief Report, we extend this scheme and in- 
troduce a general route to reduce inhomogeneous broad- 
ening effects in nanodevices. The strategy we propose 
exploits tunability to counteract the effect of noise chan- 
nels opened up by control parameters themselves. Beside 
the interest for quantum gates engineering, on a more 
fundamental level our analysis is relevant for optimizing 
fault-tolerant architectures^, showing that the influence 



of 1/f fluctuations in the solid state can be limited by 
exploiting the band structure of coupled nanodevices. 

We apply this method to a universal two-qubit gate in- 
volving an entanglement operation on two quantum bits, 
a necessary step toward the construction of a scalable 
quantum computer—. We consider a yi — SWAP gate 
based on a fixed coupling scheme and show that, for se- 
lected optimal couplings, it can be accurately operated 
even in the presence of 1/f and high-frequency noise. As 
a relevant example, we analyze a realistic charge-phase 
two-qubit gate^. 

Optimal tuning - We denote the device Hamiltonian 
with "Ho, and its eigenvalues as w,-. Fluctuations of 
the control parameters due to external environments are 
modeled by T-L\ = — \ J2 a Qa &> X a . Here {X a } are col- 
lective environmental variables coupled to the nanode- 
vice operators, Q a . Typically, power spectra of X a dis- 
play a 1// low-frequency behavior up to some cut-off, 
followed by a white or ohmic fiank 2 -^. Defocusing is 
due to low-frequency fluctuations whereas noise at high- 
frequencies (of the order of the eigen-splittings) is respon- 
sible for system-bath inelastic energy exchanges. Act- 
ing on different time scales, these two processes can be 
considered separately. Upon extending the multi-stage 
elimination approach^, we separate the effects of low- 
and high-frequency components of the noise by putting, 
X a X a (t) + X*. Here {X a (t)} = X(t) are classical 
stochastic variables describing low-frequency (1//) noise, 
and can be treated in the adiabatic approximation (adi- 
abatic noise). Instead, high-frequency fluctuations X' a 
are modeled by a Markovian bath and mainly determine 
spontaneous decay (quantum noise) . Nanodevice popula- 
tions relax due to quantum noise (Ti-type times), which 
also leads to secular dephasing (T 2 * = 2Xi-typc). De- 
focusing due to low-frequency noise determines further 
decay of the coherences. 

In the adiabatic and longitudinal approximation^ the 
system evolution is related to instantaneous eigenval- 
ues, o>i(X(£)), which depend on the noise realization. 
The leading effect of low-frequency fluctuations in re- 
peated measurements is given within the Static Path Ap- 
proximation (SPA). It amounts to replace, in the path- 
integral average over X(i) realizations, X a (t) with sta- 
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tistically distributed values X a (0) = X a at each repe- 
tition of the measurement protocol. As a consequence, 
level splittings Wy-(X) arc random variables, with stan- 
dard deviation S 



= \/<<H) - (<Jw -) s 



where <5c<j.jj = 

Wjj(X) — Wjj. The splittings enter the evolution of 
the reduced density matrix coherences in the eigenba- 
sis of Ho as Pij(t) w Py(0) / dXP(X) exp[— iwjj(X)t]. 
The probability density can be taken of Gaussian form 
P{X a ) = cxp[-X2/2E| :a ]/v / 27rSA' Q in relevant cases^. 
Optimal tuning of control parameters is achieved when 
the variance of the relevant instantaneous splittings, S^-, 
is minimized. It results in the reduction of Pij(t) decay 
due to defocusing processes. 

If Wjj(X) is monotonic in a region \X a \ < 3Ex Q , 
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The variance attains a 



minimum for vanishing differential dispersion. This is 
the case of single-qubit optimal points. For the charge- 
phase two-port architecture, control is via gate voltage, 
q x = C g Vg/{2e), and magnetic-flux dependent phase <5, 
entering the Joscphson energy Ej = Ej cos S. Thus 
X a —> Aq x , AEj and the optimal point, q x ~ 1/2,5 = 0, 
is at the a saddle point of the energy bands^. When 
bands are non-monotonic in the control parameters, min- 
imization of defocusing necessarily requires their tuning 
to values depending on the noise variances. We illustrate 
potentialities of this result for a universal two-qubit gate. 

yi — SWAP gate - A simple way to realize two-qubit 
entanglement is via a fixed, capacitive or inductive, cou- 
pling scheme^. Fast two-qubit operations and coupling 
switching on/off are achieved by individual-qubit control. 
Fixed coupling has been used to demonstrate two-qubit 
logic gates H e )ii£ and high-fidelity Bell states generation 
in capacitive coupled phase qubitsil. 

The interaction is effectively switched on by tuning the 
single-qubit energy spacing to mutual resonance. The 
building-block Hamiltonian reads Ho = — § <g> — 
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al-'csaf (h=l, <Tf'\±) =±|±». The 
Hilbcrt space factorizes in a " SWAP-subspace" , spanned 
by the eigenstates of H {|1} = ^- (-| + -) + | - +)) 

and |2) = i(| + -) + | - +))} (eigenvalues ui = 

=po; c /2), and in a "Z-subspace" generated by the eigen- 
states {|0) = — sin^ | + +) + cos^ | ) and |3) = 

cos ^ I + +) + sin ^ I )}, where tan 95 = — w c /(2f2) 

(eigenvalues wo 
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Wc focus on the 



a/i - SWAP operation | + -) -> \ip e ) = [| + -) - i\ - 
+)]/\/2 which generates by free evolution an entangled 
state at ig = ir/2u> c . We consider the general case where 

H\ = -\ xia[ 1] + ZiaP ®l 2 -|li® &2(xf ) + ha£ } . 

Since transverse (£{) and longitudinal (z,) fluctuations 
usually have different physical origin, we assume they 
are independent^. For capacitive coupled charge-phase 
qubits, polarization fluctuations cause transverse noise 
Xi oc 4E c ,iAq x ,i (Ec,i qubit i charging energy), phase 
fluctuations lead to longitudinal noise z% cx AEj^. 
Both low-frequency fluctuations induce a stochastic ef- 
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FIG. 1: (Color online) Dispersion of the SWAP and Z split- 
tings for uj c /il — 0.01. Left: Jo^i/fi from numerical di- 
agonalization of Ha + Hi (top). Zoom around the origin 
highlights the interplay of 2 and 4 order terms, barrier 
height oc cjj? (bottom). Right: Comparative behavior of dis- 
persions in the two subspaces. Top: SWAP exact splitting 
blue (gray) line, expansion |T} for X2 = 0, yi — dashed blue 
(gray), 2 nd order expansion (dotted), Z splitting (thick gray) 



5uj 30 ~ -2*2. |( yi +y2 ) + [l + I (<£ ) 2 ] }, and single 

qubit dispersion (circles). Bottom: Longitudinal dispersions. 
The Z-subspace (light gray) is much more sensitive both to 
transverse and longitudinal variations. 



fective SWAP splitting, W2i({xi, z{\), the relevant scale 
for the gate operation. It can be obtained from Ha + Hi 
treating perturbatively the stochastic fields in Hi, 



V2i(xi,x 2 ,zi,z 2 ) ~ uj c - 777^(^1 + -J— ( z i - z a) 2 

+ ^j5(*? + 4)(*i + «) + - *a)(*i - a) (1) 



A key feature is that lo 2 i is non-monotonic in the small 
coupling uj c <C SI. This is due to a selection rule for lon- 
gitudinal fluctuations. They only mix states inside each - 
SWAP or Z - subspacc, while a^-fluctuations mix the two 
subspaces. For instance, second order transverse correc- 
tions to uj\ are 2 I (-M^iN) | 2 /( w i ~ UJ i) °t w ci whereas 
longitudinal ones vary as w" 1 oc |(l|'Hi|2)| 2 /(wi — W2). 
Non-monotonicity in uj c results in a competition be- 
tween 2 nd and 4 th order Xi-terms in Eq.fT]) and in non- 
monotonic band structure, Figure 1 (left). Because of 
this subtle feature, identification of the best operating 
condition necessarily requires consideration of the noise 
characteristics. Indeed an optimal coupling can be found 
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FIG. 2: (Color online) Envelope of the concurrence in the 
SPA, SI = 10 1 rad/s. To elucidate the significance of the 
optimal coupling scheme, we consider 1/f amplitudes larger 
than expected from single qubit measurements^. Left: Ef- 
fect of transverse noise with T, x /Q, = 0.08 and uj c /CI — 
0.01,0.02,0.06,0.08,0.1 (dotted, dot-dashed, triangles, thick 
(gray) red, dashed). Light gray is the single qubit coherence, 
\p+-(t)\ = [1 + (£^/n 2 )]- 1/4 £. Inset: C(t) and its envelope 
for coc/n = 0.02, 0.08. For optimal coupling u c ~ Ex, at 3ns 
already 8-SWAP cycles occurred. Right: Effect of transverse 
plus longitudinal noise on qubit 2, T, Z2 /Q = 2.5 x 10~ 3 . Inset: 
Effect of longitudinal noise, u c values as in left panel. 



which minimizes the SWAP variance 




(2) 

where we assumed equal transverse variances, T, Xi = T, x 
and T, Z1 <C S 22 , which mimic typical experimental con- 
ditions^. Eqs. (HJ) - ([2]) highlight the general result of 
this Brief Report. While higher stability with respect 
to longitudinal fluctuations is attained by larger cou- 
plings, minimization of the detrimental transverse low- 
frequency noise components is obtained by tuning the 
coupling to an optimal value Q c . For S Z2 <C S x , this 
is the transverse noise variance, Ca c w 2 1 / 4 Y, X . It can 
be estimated by independent measurement of the ampli- 
tude of the 1/f transverse noise on the uncoupled qubits, 

sl /f = 7r£2[i n(W7m)wr i ; E 2 = r-^/^H 

(low and high frequency cut-offs 7 TO and 7a/). Note the 
higher stability of the SWAP-splitting compared to the 
qubits Larmor frequency and the Z-splitting, Figure 1. 

Working at the optimal coupling minimizes defocusing 
and guarantees excellent performance of the Vi - SWAP 
operation. As a unambiguous test of entanglement gen- 
eration and its degradation due to noise, we calculated 
the concurrence during the gate operation^. If the sys- 
tem is initialized in | H — ) , the 4x4 density matrix in the 
computational basis is non-vanishing only along the di- 
agonal and anti-diagonal at any time (X-states) and the 
concurrence takes a simple formal. 

We first consider the effect of low-frequency fluctua- 
tions. In the adiabatic approximation, the concurrence 
simplifies to C(t) = 2\lm{p 12 (t)}\. With the SWAP- 
splitting expansion ([T]) (including in the 3 rd and 4 th order 




V 20 40 60 80 100 20 40 60 80 



t (ns) t (ns) 

FIG. 3: (Color online) Concurrence for ui c /Q. = 0.01 blue 
(gray) line and for optimal coupling Cu c /Q. — 0.08 (light gray). 

(a) Effect of high-frequency noise, S Xi (<*>) ~ 8 x 10 s s _1 , 
S Z2 {uj) « 4 x 10 7 s" 1 > S z1 (oj c = 0.08f2) « 6 x 10 3 s" 1 
— . Inset: At short times C(t) w 2|Im{pi 2 (£)}| (diamonds). 

(b) Effect of 1/f noise (parameters as in Fig. [2| and white 
quantum noise. Inset: Asymptotic behavior. Results are 
minimally modified considering the dynamics of fluctuators 
generating 1/f transverse (longitudinal) noise in y m = 1 s _1 , 
jm = 10 6 (10 10 ) s -1 (numerical solution of the stochastic 
Schrodinger equation). 

the terms cx ui^ 1 ) in the SPA we obtain 

Pi2(t) = P u(0) ^y^e-W) K [h(t)] (3) 

where h(t) = (E^ + 53^ + iuJc / t ) (ft 2 /T? x + iu c t) 2 / '(4ft 2 ), 
and -Ko[^] is the K-Bessel function of order zero—. By 
increasing the coupling to match the optimal condi- 
tion two goals are simultaneously achieved: minimiza- 
tion of initial defocusing and fast two-qubit gate (Fig- 
ure [2]) . The first SWAP error takes remarkable values 
e = l-{i> e \p{t E )\i>e) « lCT 3 -lCr 4 , for uj c a Z x < 0.05ft 
(numerical simulations). This is an interesting effect, 
considering that single qubit coherence times would be 
rather small at the same 1/f amplitudes, T 2 ~ 5 ns. 
The optimal coupling scheme is effective against large 
amplitude 1/f - noise even in the presence of high- 
frequency fluctuations. Within the secular approxima- 
tion, quantum noise leads to additional exponential de- 
cay of SWAP-coherence, p\ 2 {t) = p\ 2 {t) exp{— T 12 t}. 
Thermal relaxation processes (fc^T <C ft) populate 
levels i = — 2 and the concurrence reads C(i) pa 
V(Pii-P22) 2 +2(Im[p 12 ]) 2 - I sin^poo. The SWAP de- 
cay rate is related to escape rates from levels 1 and 2, 
f 12 w where T\ = T w + T 12 , T c 2 = T 2Q + T 21 and 

T l0 cx [S xl (oj i0 ) + S X2 (u i0 )], r 2 i cx [S z1 (oj 2 i) + S Z2 (u 2 i)]. 
These rates enter the populations in the combinations 

r± = -(r? + r|)/2 ± v /(r5-r°) 2 + 4r 12 r 21 /2S. The 

SWAP-coherence rules the relevant short-time behavior, 
jpia| oc 1 - T 12 t (or |p ia | cx 1 - {T, 21 t) 2 /2) depend- 
ing on the most relevant quantum (or adiabatic) noise 
component. Exponential long-time decay is due to the 
strongest longitudinal or transverse quantum noise, Fig- 
ure [3J The finite asymptotic value reflects the entan- 
gled thermalized state (no "entanglement sudden death" 
occurs^). The above entanglement characterization 
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FIG. 4: (Color online) Qubit 1 switching probability to |— }, 
P m (t) and probability P (2) (t) to find qubit 2 in the initial 
state |— ) in the presence of 1// and white noise for lj c /Q = 
0.01 blue (gray) line and for optimal coupling u c = — 
0.08O. (light gray), (a) P^ 1 ': Exponential short-time limit at 
u) c , algebraic otherwise, (b) P (1) and P (2) [dashed (gray) red] 
anti-phase oscillations for w c (main), u c /Q. = 0.01 (inset). 

translates into directly measurable quantities. Out of 
phase oscillations of single qubit switching probabili- 
ties signals two-qubit states anti-correlations 1 ^. It fol- 
lows from P(a)(t) = p(t) ± Re[pi 2 (i)], with P(t) = 
— \ cos^ [pn(i) + P2i{£)\ +cos 2 (y). For a charge-phase 
y/i — SWAP gatei^, defocusing due to 1// polarization 
and phase noise is considerably reduced at optimal cou- 
pling, Figure 2] Phase quantum noise on qubit 2 dis- 
placed by its optimal point contributes to initial decay. 
Oscillations visibility is larger than 90% until w 10 ns, 
corresponding to w 25 SWAP cycles. This contrasts with 
strong initial decay for non-optimal coupling. Long-time 
exponential decay is due to polarization quantum noise. 

In conclusion, we have proposed a general route to re- 
duce inhomogeneities due to 1 // noise by exploiting tun- 



ability of nanodcvices. As a relevant and timely issue, 
we have illustrated the considerable improvement of the 
entangled dynamics in a universal two-qubit gate. It re- 
quires preliminary noise characterization and work in the 
protected SWAP-subspace. In the considered scheme, 
coupling is controllable by dynamically tuning the qubit 
frequencies. We demonstrated that an efficient gate may 
be obtained with no additional dynamical decoupling 
protocols, even if at least one qubit has to be moved away 
from its own optimal point. In order to implement a scal- 
able architecture, this scheme has to be supplemented by 
switchable coupling between small sub-units, formed by 
single qubit and universal two-qubit gatesi^. Note that 
our approach can also be applied to directly switchable 
coupling schemes which have been recently proposed as 
alternative, potentially scalable, design a 24 ' 25 . 

In the broader perspective of fault tolerant architec- 
tures in the solid state, our work provides a strong hint 
on how to extend the analysis of Ref.— to include the 
effect of time-correlated phase errors, typically affecting 
nanodcvices. Note, in this respect, that fixed coupling 
schemes are at the basis of qubit encoding^ 1 -. 

Finally, we remark that the device reliability may be 
qualified by the impact of one/few impurities strongly 
coupled to These induce non-Gaussian fluctuations 

which randomly displace qubits from resonant condition, 
possibly resulting in limited readout fidelity^. These 
effects are neglected the present analysis, but could be 
addressed as suggested in Ref4. 
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